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CURRENT TO MOVING SPHERICAL AND CYLINDRICAL 
ELECTROSTATIC PROBES 

ABSTRACT 

The current collection characteristics of moving cylindrical 
and spherical electrostatic probes are evaluated for velocitir , en- 
countered in earth satellites and planetary probes. Simple asymp- 
totic current formulas are developed for large values of the Mach 
number and potential. For Mach numbers greatei than 2.5, the or- 
bital motion limited current is given within 1% by 

d - 1 

i + JlL + 2L\ 2 

mw 2 mvv 2 / 

where d = 1, 2, 3 for a plane, cylinder, or sphere respectively. A 
data reduction technique for determining ion densities from the ac- 
celerated ion current is discussed. The current expressions are 
given in terms of standard functions and recurrence relations, 
therefore allowing straightforward computer programming of all 
the probe current formulas. 
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INTRODUCTION 


Cylindrical and spherical electrostatic probes have lx»en widely employed 
on rockets and satellites for charged particle temperature and density measure- 
ments. 1 ' 3 In the future these devices may be used on planetary missions. A 
common feature of the planetary and earth electrostatic probes is the rapid 
vehicle velocity relative to the plasma thermal velocity. This results in high 
ion Mach numbers and non-negligible electron Mach numbers. 

The formula for the current to a moving cylinder has been obtained bv 
Kanal 4 ; however the series expansions derived in reference 4 are not valid for 
high Mach numbers. A number of authors s have investigated the current to 
a moving sphere, however only Kanal has obtained the current expression for 
an arbitrary sized sheath. 

We reexamine the expression for the current to a moving cylindrical probe 
(eqs. 15 and 27 of reference 4) in order to derive expansions valid at high Mach 
numbers. YVe include the expressions for the current to a moving sphere so 
that the asymptotic current expressions may be compared with those of the 
cylinder. 

The key assumptions used in deriving the probe currents are that: 1) the 
probes are surrounded by a sheath of ideal geometry (cylinder or sphere); 2) 
no collisions take place within the sheath; and 3) the particle distribution out- 
side the sheath is Maxwellian with a superimposed drift velocity. The sheath 
symmetry becomes distorted at high Mach numbers; however most of the col- 
lected current then originates in a small region about the velocity vector where 
sheath symmetry is approximately maintained. The effects due to the finite 
length of the cylinder probe will be neglected in the present treatment. 
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BASIC CURRENT EQUATIONS 


The derivation from first principles of the integral expression for the 
current to a moving sphere or cylinder has been treated in the literature. Both 
the sphere and cylinder probes may be treated by the method of Kanal. 4 We ex- 
press the current formulas in normalized form, 


I = i/i 


random * 


( 1 ) 


where i is the actual current and i raM(|t)in is the random current to an uncharged 
probe of area A due to particles of mass m with number density N and at temp- 
erature T, 


i 


random 


■( 


k T 
2 77 m 


1 / 2 


NAe 


( 2 ) 


The normalized accelerated particle current is given by the formula, 


d - 1 


sds (s 2 + V) 2 It'D ( S ,M) 


/v 


t 2 


(?) 


d - 1 


/•yv 


1 / 2 


sds s d * 1 I (d * (s, M) , 


(3) 


where d = 2, 3 for cylinder and sphere respectively, and 


2 2 


1(2) (s, M) =-l— e - (M ‘ +s >I 0 (2Ms), 


77 


1 2 


(4) 


I (3) (s, M) = - 


0-(s -M) 2 e - ( s + M ) 2 


2 M S 


( 5 ) 



The symlxds are defined by: 


M (m w 2 sin 2 f ; 2kT)* 2 Mach number, 


0 angle between probe axis and velocity 
vector for cylinder; 0 v 2 for sphere, 

w probe speed relative to the stationary plasma, 


y = l/( t 2 r 2 - l) 1 2 , 


a shea t h radius , 


r probe r adius , 


V 


e P 

kT 


ratio of potential to thermal energy. 


The normalized retarded particle current is: 

od d - 1 

I rt>t r I sds (s 2 - V) 2 I (d) (s, M). (t>) 

J v > 2 

In formula (4), I 0 is the modified Bessel function of order zero. 8 The in- 
tegrals for the sphere may be readily performed, and yield exponentials and 
error functions. The integrals for the cylinder can at best only be reduced to 
the form of a power series whose coefficients are known functions. Kanal 4 ob- 
tained a single power series in powers of M/V 1 2 with coefficients depending on 
Bessel functions and incomplete gamma functions for I and modified Bessel 
functions for I ^ . 
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We derive a power series representation for the cylinder current in power 
of M 2 and V from which asymptotic formulas valid for high M or V can lx» ob- 
tained by summing over M or V. The summation yields series in confluent 
hypergeo metric functions. 8 

The evaluation of the cylinder current reduces to the problem of performing 
the integral, 


I (x. y) 


rr 


4_ 
1 2 


L s 


sds (s 2 f y) 1 2 ♦ M ‘> I 0 (2Ms) 


( 7 ) 


The dependence of I (x, y) on M 2 is simplified by performing a lAplace transform 
in M 2 : 


i I (x. y) 



d M 2 e“ pM I (x, y ) 


<«) 


Using the Laplace transform of the modified Bessel function \ we obtain the 
result. 


Cl(x.y) 
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1 2 1 + p 


ds(s + y) 1 2 e ~ ps + p) . 


O) 




This is in the form of an integral representation of the confluent hvpergeometric 
function 8 , / (a, c; z). Therefore using the relation, 

yp (a, c; z) z 1 " c 0 ( a-c + 1, 2-c;z), 0^) 


t 

t 


/ 


4 


4 


we find: 


I (x. y) 


~ 1 • p* ( 1 + p ) / 1 * P \ 

trl 2 Up \ P / 


/ I ! p(x • y >\ 

V 5 ” 5 ’ l,p ) 


01 ) 


To facilitate the separation of variables and taking the inverse Laplace 
transform, we partition I into two parts, 

I (x. y) I, (y) ^ I 2 (x. y) 


(12) 


I (x. y) = £ I, (y) ♦ £ I 2 (x. y), 


according to the relation between , and i , 


0 (a. c; z) r , ( 1 ~ i (a, c; z) 

I (a - c i 1 ) 


— — z 1 ~ c (a - c + 1. 2 - c; z) , 

r (a) 


( 13 ) 




T 

a, c; z) = y. 

II r 0 


( C )„ 


(!-») 


The result of this partition is: 




1 /I +P 


3 2 


»py (1 + p ) 


P V P 


( 15 ) 


£ I 2 (x. y) = - f-L (* + y> 3/J e- 0 (l. |: £^f)/r (f) . 


( 16 ) 
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In performing the inverse Laplace transform rn Equations (15) and (Hi), the 
path of integration reduces to a circle of radius greater than unity and centered 
at p = -1. For I, (y) a branch cut extends from p = -1 to p = 0. Both integrals, 
1, and I , are evaluated with the change of variable p ■ -1 ♦ q, and an expansion 
of the integrands in powers of q. Evaluating the residue at q * 0 we find: 


I, (y) 


L L 


( - M 2 ) n y m r (n ^ m - 1 2) 
n! m! I~ (m - 1 2) P (n * 1 ) ' 


07) 


I 2 (x. y) 


-(x + y) 3 2 




(-M J ) n (-X)" (x + y) r (n » r i 1) 
n! m! ♦ -5. \ (m ♦ r ♦ 1 ) i (n * 1) 

(18a) 


I 2 (x. y) = - (x ♦ y) 3 2 


X <T 

L E 

m s 0 n a 0 


(-x) m (x + y ) n / (lunul, 



1; -M 2 ) 


(lsli) 


These basic formal power series will lx? rearranged later as series in 
confluent hypergeometric functions (in the manner of Equation 18b) which can 
either be programmed on a computer o expressed in analytic, asymptotic form. 


Finally we give the connection between the basic formulas, Equations l 7 
and 18, and the normalized cylinder currents. 


1.CC = 1.<V> * I 2 (/ 2 V, V) -2.1, (v* V. 0) 

c y I i nHe r 

*,.« =1. (-V). 

cy l i nde r 


(19) 

(20) 


We note that for sufficiently small values of M and V, the power series 
representation of Equations 17 and 18 are rapidly convergent and would there- 
fore be appropriate for computer evaluation. 
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Two limiting forms of the accelerated current are commonly employed: 

1) the sheath area limited case when the potential is small and the sheath collapses 
to the prolx?, a/r • 1 and / • ' ; and 2) th-, orbital motion limited case when the 
sheath is much larger than the probe, a/r 1, > • r/a. The normalized current 
expressions for these two ideal cases are: 



cy I i ndp r 

[ o ml Ij (V) + Ij (0, V). 

c y I i rule r 

The expressions for the normalized current to a sphere are obtained by 
integrating Equations 3 and (i: 


(21) 

( 22 ) 


^IICC 
s phe re 


2 ( M '4) 


2 erf (M) 


I ) 

+ — [ e r f (M - y V 1 2 ) + e r f (M + y V 1 2 ) ] + -L e ' M ^ 




(23) 


r 


- 1 < 


, 1/2 


M 2 


H 


2 M 


[erf (M - y V 1 2 ) + erf (M + yV l : )] 






M + yV 1 2 - ( m - > v* V , M - /V 


4 M 


. 1 2 »-(M + rv' 2 ) 2 

4 M 


> . 




I 
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1 2 


M 2 


H 


ret 2 

s phe re 


2 M 


[erf <M-V l 2 ) + erf (" + V 1 2 )1 


7 





( 24 ) 


M ♦ v * 2 v -(M .V 1 7 ) 2 

4M 

, M - Vl 2 ^(M + v’ V 

4M 

For the two idealized cases the accelerated current reduces to: 


(25) 

(26) 

Equation 25 has been given by Kanal 7 and Equation 25 by Sagalyn. 6 


I 


a 2 I - 1 2 


sn I 


Hphero 


M J ♦ i "] 

■jjj— erf (M) . 


7 T 


1/2 M 2 


OKI | 2 

s phr r p 


M 


erf (M) fi-e' M 
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SHEATH AREA LIMITED CURRENT 

The sheath area limited current is determined by the condition that all 
particles that enter the sheath are collected; it is obtained from the general 
accelerated current by taking the limit as a/r * 1 and multiplying the result by 
( a/r) d ~ l . 


For the cylinder, the sheath area limited current is given in closed form: 

I.., ^ (-1. 1; . (27) 

cy 1 i nde r 

Using the contiguous functions, the expression for the derivative of t and 
the confluent hypergeometric representation of the modified Bessel function, 
one can demonstrate the equivalence of this result with Equation 22 of Reference 
4. 
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I 


Sll I 


cy I » rule r 


L<-- sr/i a 


• M 1 ) I 0 






A comparison of these two representations demonstrates the usefulness of 
the confluent hypergeometric functions; only one function need be evaluated in- 
stead of an exponential and two Bessel functions and the asymptotic form for 
large M is easily obtained from /. For small M, the Taylor series representa- 
tion of t provides a convergent calculation of I t . 


For large M, the asymptotic series is derived from the following exact 
series which is obtained from the integral representation of /: 


t (a. c; - 7 ) 


uu 

L 


n - 0 


(n - c + l) n I (c) > (a + n. z) 
n! I (c ••) (*) *• ♦** 


(28a) 


a - c i nega t ive in t e ge r , 


t (a, c; z ) 


‘ E 


II ; 0 


(1 -a) n 1(c) y (c - a +n, z) e* 
n! I (a ) I *(c - a ) z r " • + n 


(28b) 


a / negative integer, 


where the incomplete gamma function may be written as 

> (a, x) - a ' 1 x a (a, a + 1; -x) 


For sufficiently large z, the following equation holds; 

1(a) - y(a,z)' v 'z n ’ 1 e~ z , 


Therefore we have the asymptotic expansion: 


9 




(M InrRc ) 


Writing out the first few terms we have, 


s a I 

cy I i nde r 


n 

r 


1 


1 


,1 /2 


V 

( M large) 


4 M 2 32 M 4 



(30) 


Retaining only the iirst term and multiplying by the random current, we 
find that the actual current is equal to 



sin 6 

V 


» 


which is the product of the number density, charge, relative velocity, and the 
cylinder probe area projected normal to the velocity, 


-lL A sin 6 \/v - 2 a i sin Q i . 
r 

The bphere sheath area limited current is given by Equation 25, the asymp- 
totic form for large Mach number, M, is: 


sa 1 

sphere 


^ 77 


1 / 2 


M 


(M large) 



(31) 


where we have neglected terms of order e 
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RETARDED CURRENT 


The retarded current is independent of the sheath size and depends only on 
the potential across the sheath. The cylinder retarded current is given bv the 
double series (Equation 20): 


I 


ret 


c y 1 mder 


« 00 

L L 

i) j . 0 m - 0 


(- m 2 ) n (- vr r 

( 1 ^ 

\ 2 j 

1 

n! m! r ( m - j) 

: (n + I) 


(32) 


For small M and V this series provides a convergent means of calculating the 
current. 


For small M and arbitrary V, we sum the series in V giving: 


cy i i ndc r 




(-M 2 )" 





(33) 


The coefficients of the powers of M 2 {the t) are polynomials in V multiplied by 
e v . Knowing the first two functions: 



4 > Q. - = (1 + 2 V) e- v . 

the general terms can be found from the recurrence relation: 



(34) 


( 35 ) 
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This provides a useful computer scheme for evaluating Equation 33. Writing 
out the first few terms, we have 


cy 1 i ndo r 

The series given by Equations 33 and 36 would adequately treat the retarded 
electron current for earth or planetary probes. 

For higher Mach numbers we perform the sum over M 2 to obtain: 


1 + 2 V 


M 2 - 


(1 + 4 V - 4 V 2 ) 
16 


M 


-} 


( 36 ) 


I 


... 

cy 1 inde r „ * o 


(37) 


The first term in this sum is proportional to the sheath area limited current; 
the general terms are determined from the first two 4 and the recurrence 
relation: 


4> f-i. l; -M 2 


) =e- M ' 2 (1 4 M 2 ) I 0 4 M 2 I, 

— 

1; -M^) - e-' 2 I 0 (*) 

(m 4 1; - M 2 ^ = U2m - M 2 ) 4 > 4 i, 1; - M 2 ^ 


M 2 

2 



(38) 


Using the asymptotic form of 4 for large M, this series becomes, 


05 00 


cy I i nde r 


4 L L 

n = 0 m = 0 


-V \ m 

M 7 J 


M 


LV- (..* 


' 1 

m 

2 


n! m! i - m 


(39) 
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The summation over the index m converges for M"’ ' V and gives a series in 
1/M 2 : 




cy l i nde r n * 0 


C 0 — (M J - V) 1 
'J 


C . _i M J _ C . 

"♦! dM a " 


( 40 ) 


The first few terms of the asymptotic expansion are given as: 


I 


ret 


c y I i nde r 


,1 2 


M 


M > V 


',-^V 

, M J I 


l 

2 M 2 


\ X. \ ‘ l 2 1 (\ IV 3 2 

M 2 ] 2 l M 2 / 


1 

4 M 2 


i.iV".3 A.x 


5 2 


M 2 


M 2 


(41) 




» M 2 / 


The retarded current to a sphere is given by Equation 24; the asymptotic 
form for lar ge M is: 


I 


77 


1 2 


ret 9 

4 m 

sphere M large 


M 1 + 


1 V_ 
2 M 2 M 2 


(42) 


ORBITAL MOTION LIMITED CURRENT 

The orbital motion limited current is independent of the sheath size and 
depends only on the potential across the sheath. The current to a cylinder 
(Equation 22) is given by the double sum: 


13 




r 


no no 


I 


om I 



(-M 2 )" V m 


cylinder „ = 0 m 0 


n! m! 


< 


v. 


1 ( nf "-j) 

\ " 5 ) I ' ( n ♦ 1 ) 1 + 4 ) 1 < n ♦ ! 1 > 


V 3 2 I (lumtl) 


> (43) 


The formula suitable for low Mach numbers is obtained by summing over V: 


1 . =— - 

° ml I 2 

cylinder 11 n „ 0 


E O.V'K-H 


(44) 


From the recursion relations among the ^ and the fact that, 

0 - 77 l 2 e v erfc V* 2 , 

/% 

we can determine all the coefficients of M : 


0 


(-1, -I; ILL! c v erfc V ,/2 + V 1 2 , 


(45) 


(46) 




e v (1 - 2V) erfc V 1 2 + V 1 2 , 


0 


-i, *) .[(2„.|.V)* -1;,) 


-0 


( 1 1 

\ 2 2 



(47) 


(48) 


(n + 2) 


Use of the recurrence relation (48) in Equation 44 provides a convenient 
means of evaluating the orbital motion limited current to a cylinder for M < 4. 
Writing out the first few terms we find: 


I . = ^ v-V + e v erfc (' V) + — 

om 1 — ' ' 2 

cylinder 1 ,r 

M smn I I 


-L^V + (1 - 2 V) e v erfc (vV) 

V 77 


M 4 

8" 


2 

(V + - Vv + ( 

_ V v 

\ 2 ) { 


1 - 2 V - 2 V 2 

2 


) 


e* v erfc (v V) 


(49) 
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Using the exact representation for , : 


'/'(a, c: 7 ) / 

t—i n! (n) 

n r 0 


/ (a + n. z ) 
(-z) n z a 


1 - c 


(- z ) n (c - 1 - n. ?. ) 


(50) 


we obtain an asymptotic expansion for . at large potentials: 





2 -n 



(f i) 


Therefore we obtain an expression valid lor V » 1: 



cy I i nHe r 


t l 

m s 0 n r 0 





When V > M J the series in n converges and is given by: 


om 1 

c y I l rule r 



-4: /V + M 2 . 

v v 


n + 1 



M 2 3 \ dc n 
ntl 3M 2 / v 


(52a) 


A. I 

om 1 

c y 1 inder 


^ 2 
_r_ v 1 / 2 

V 77 

V > M* 





1 

4 V 2 


\ ' 3 ^ 9M^ L M 2 \ 

V / ' 4V \ V / 


(52b) 


15 M 4 

8 V 2 



15 


Wc obtain an expression useful for high Mach numbers by summing over M : 



The second set of confluent hvpergeometric functions, t ( m ♦ 1, 1; - M 2 ), 

2 

consist of Lagv.orre polynomials in M 2 multiplied by e ’ M , therefore we may 
neglect them in the asymptotic expansion. Consequently the orbital motion 
limited current for large M is identical with the asymptotic retarded cylinder 
current with the substitution V - - V: 


I 


om I 


cylinder M I n r ' 




1 2 


1 

2 M 2 


i M 2 






It is interesting to note that up to the second term the expansion for I . , 

* r om I ’ 

when M 2 V and M 2 > V, Equations 52 and 53, are identical. 


Taking the lowest order terms we find that Equation 54 can be approximated 
by the single square root expression: 


I 


om 
c y I 


2 / l 

I " — 1*2 

> nde r M‘ > V 


M 2 -r V 


1 2 


(55) 


In Figure 1, there are displayed contours of the relative error between 
approximation 55 and I on , computed either from Equation 44 (for small M) or 
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Equation 54 (for large M). It is seen that the relative error is less than 1% for 
M s 2.5 and all values of the potential. 


From Equation 52, the asymptotic expansion for M‘ 
approximation: 

t.„. *-L<i + v ♦*)*'* 

cy I i nde r 1 rr 


\\ we ol)w.iin the 


(5<i) 


In Figure 2, there are displayed contours of the percent relative difference 
between approximation (58) and I , as computed by Equations 44 and 54. For 
potentials V s 5.5 and any value of M, the relative error is less than 1%. We 
see from Figures 1 and 2 that approximation (55) is most appropriate for fixed 
M when the potential is varied, and that approximation (56) is most appropriate 
for fixed V when M is varied. It is pointed out that Equations 55 and 50 are the 
generalizations to lx» taken for the Mott-Smith and Langmuir small cylinder 
current expressions 10 : 


om 1 

cy I » nde r M = 0 


_i v 1 * v . 

V7T 


Mott-Smith and Langmuir 


I = 

om 1 

cylinder M > > 1 


_L ^vTP , 

Vv 


Mott-Smith and Langmuir 


(57) 


(58) 


Where Equation 56 is the generalization of Equation (57) and Equation (55) 
is the generalization of Equation (58). It is interesting that the thermal effects 
are added to the Mott-Smith and Langmuir formula (58) by adding the term, 
1/2, under the square root. 


The orbital motion limited current to a sphere is given by Equation (26); the 
asymptotic form is simply: 
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orn I 

(iherr V 



1 

2 M 2 



neglec ting terms of order e 


( 59 ) 


Thus we lint! that the approximation for large M 2 for the actual orbital 
motion limited current to a sphere and cylinder (and plane) are similar in form, 


I , = New A 

on 1 

II *> 1 


kT 

m w 2 


ci -1 



( 50 ) 


Where A is the probe area projected perpendicular to the velocity vector and 
d = 1,2, .’1 for a plane, cylinder, and sphere respectively. 


ACCELERATED CURRENT 

The accelerated cylinder current as given by Equation U> is the sum of 
three terms, one double sum and two triple sums. We will not write down the 
formula as it is contained in Equations 17-1‘J. For small values of V and M 
this power series would be suitable for computer evaluation. For small M and 
arbitrary V we may express the accelerated current as a power series in M 2 : 


I 

acc 

c y I i nde r 



oc 


L 

n - 0 


(M 2 )" 
(n! ) 2 




(fil) 


where the initial values of f M and g n and their recurrence relations are determined 
most readily from the integral representation, Equation (.‘1): 


f 

n + 


2 



V (n t 1 ) f n + 2 ((1 t > 2 )V) J 2 e-> 2v (, 2 V) n 


(«2) 
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f 


1 




♦ /' ) V 


l* 


•y 


2 


v 





(1,4) 


g 


n ♦ 1 



R " <r>V)"** 2 e-> av . 


g 


o 


erf (y vV) - /vV c" v 


(65) 


( 66 ) 


This scheme is recommended for computer evaluation of I f or smal 1 M, 
since it Involves only multiplication and addition after the initial terms f 0 , fj , 
g 0 have been evaluated. We write out explicitly the first two terms in the ex- 
pansion: 


I 

acc 

cy I i nde r 



M « 1 





) 



~ t -i \ 


, / \“ 

N 

M 2 

U-\) 

e v e r f c /(I f > 2 ) V 

f — — e r f ( /y 2 V ) 

. . . ► 


L\-> / 


2 r v ' 

J 


For large M, Equation 18b) together with the asymptotic form Equation 28b), 
show that 1 2 -0. Therefore the accelerated cylindrical current for large M is 
I j (v) which is given by Equation (54) or (38) with V- -V. 

I = I, (V) = I . for M 2 » 1 (68) 

acc 1 v ' oml 

cy 1 i nde r 


The accelerated current for a sphere is given by Equation (23); ior large M we 
have, 


* * _ i 

acc oml 

sphere sphere 



t 


1 

2 M 2 



M 2 » 1 


( 69 ) 
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Therefore it is a property of probes with simple geometry that the accelerated 
current is given by the orbital motion limited current at large Mach numliers. 


ION DENSITY DETERMINATION 

We consider a methcxl of determining the ion densities from the volt-ampere 
characteristics of a moving probe. We treat cily the case when the ions are 
accelerated and the electron current is negligible. Unless the probe is driven 
to very high negative potentials with respect to the plasma, it will in practice 
be difficult to obtain more than two averaged parameters from the experimental 
accelerated ion current. However let us consider the general technique of fitting 
the accelerated ion current to a polynomial in the potential regardless of the 
information content of the data. One advantage of the polynomial method is that 
the number of terms in the fit is determined only by the amount of information 
contained in the data, independent of the number of ion species present. For 
sufficiently high Mach number M 2 > V or (mw 2 sin 2 * )/2> e;, the orbital motion 
limited current is the polynomial: 


1 = 


e 2 r t w s in O' 


z 

!. 0 


a V'\ 


(70) 



m 


( s ) 


T 


N /(mV 


(72) 
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where N , is the number density and m , the mass of species a. The coefficients 
of the powers of V are themselves expansions in the weighted inverse powers of 
the mass, m s) . The highest power of V, V 1 ’, depends on the potential range 

over which the data is taken, and can be estimated by the requirement that 
/ 1/2\ (V/M 2 ) v be not much less than unity. The termination of the series for 

\P / 

a t in powers of 1/M 2 is determined by the value of the smallest Mach number. 

It is assumed that the area of the probe and w sin are known, and that T can 
be determined (possibly from the electron temperature). Then by determining 
the coefficients, a f , from a least square polynomial fit to the data one obtains 
a set of the weighted inverse powers of the mass, N, m' 1 ’ , m <2) , . . .by solving 
the set of equations, Equation 71, starting with a p . If the ion species are known 
and a sufficient number of the m s are evaluated, then the mass concentrations 
can be found. If two known species are present, then only a knowledge of the 
first two quantities a 0 and a, are necessary to determine the concentrations. 


The orbital motion limited current to a sphere at sufficiently high Mach 
number is: 


w (a 0 + a x V) , 

(73) 

t k T (1) 

v'4) 

w 2 


2 k T (1) 

(75) 

w 2 



Therefore the polynomial fit for the sphere is nearly the same as that for the 
cylinder, however with truncation of the series in V and l/M 2 . 
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An alternate method of fitting the accelerated ion current proceeds on the 
assumption that the species as well as w sin - and T arc ’.nown. Then using the 
asymptotic expansion for the orbital motion limited current, Equations 54 or 52, 
the current is a linear function of the number densities, N , and may be statisti- 
cally fitted to the data points, 


i 



f (V) N 


(7b) 


where the coefficients f (V) depend on the species (mass), potential, temperature, 
and relative velocity w. This technique works best when there are primarily 
two species present. Although it is not unique, this technique can also be used 
to guess the masses of two species by examining the sign of the number densities, 
N computed for an assumed mass pair and then comparing the net density N «* 

Nj N, with the electron density. Additional information to aid in the ion density 
determination is contained in the altitude variation of the ion species. 
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Figure 1. Contours of relative difference between the cylinder current 

approximation, — M (1 + 1/2 M 2 + V/M 2 ) 1 2 , for M 2 > V, iind the 

V T 

orbital motion limited current. 

Figure 2. Contours of relative difference between the cylinder current 

2 

approximation, -^=-(1 + V + M 2 ) 1 2 , for M 2 v ' V, and the orbital 

V V 

motion limited current. 
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gure 1. Contours of relative difference between the cylinder current approximation, 
for > V, and the orbital motion limited current. 





alive difference between the cylinder current approximation, 
motion limited current. 


